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Q> Abstract 
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J> Peak amplitude measurements of the fundamental mode of oscillation of a suspended 

IT^ aluminum alloy bar hit by an electron beam show that the amplitude is enhanced 

by a factor ^ 3.5 when the material is in the superconducting state. This result is 



consistent with the cosmic ray observations made by the resonant gravitational wave 
detector NAUTILUS, made of the same alloy, when operated in the superconducting 
state. A comparison of the experimental data with the predictions of the model 
describing the underlying physical process is also presented. 
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1 Introduction 



In a pioneering experiment [1] B.L. Baron and R. Hofstadter measured me- 
chanical oscillations in piezoelectric disks when penetrating high energy elec- 
tron beams impinged on the disks. The authors outlined the possibility that 
cosmic ray events could excite mechanical vibrations in a metallic cylinder 
at its resonant frequency and that they could represent a background for 
experiments aimed at the detection of gravitational waves (gw). The gw res- 
onant detector NAUTILUS, a massive (2.3 t) suspended cylinder made of 
an aluminum alloy (A15056) that can be cooled down to the thermodynamic 
temperature of 0.1 K, has been equipped with a cosmic ray detector to study 
the interactions due to cosmic rays and to provide a veto against the induced 
events in the antenna. The results on the cosmic ray observations made by 
NAUTILUS can be summarized as follows: 1) when the antenna was oper- 
ated at a temperature T = 0.14 K, well below the transition temperature 
from normal-conducting (n) to superconducting (s) states of the material, the 
rate of high energy signals due to cosmic ray showers was larger than the 
expectations based on the model describing the underlying physical processes 
[2,3]; 2) there was no evidence of this feature when the antenna was operated 
at T = 1.5 K, well above the transition temperature [4]. From one side the 
hypothesis that this behavior was linked to the conducting state of the an- 
tenna and on the other side the incomplete knowledge at very low temperature 
of the thermophysical and thermodynamic parameters needed by the model 
have motivated an experiment (RAP) to measure the longitudinal oscillations 
of suspended cylindrical bars exposed to electron beam pulses of controlled 
energy and intensity. The experiment, performed at the Beam Test Facility 
(BTF) [5] of the DAFNE ^-factory complex in the INFN Prascati Laboratory, 
has already obtained the following results: 1) the measurements over a wide 
temperature interval (4.5 K < T < 264 K) on a bar made of the same alu- 
minum alloy as NAUTILUS have confirmed with good precision the validity 
of the model [6]; 2) the measurements on a pure niobium bar operated in the 
n and s state have demonstrated that the oscillation amplitude of the bar 
induced by the interaction with the beam depends on the state of conduction 
of the material [7] . 



In this letter we report on the measurements made on the aluminum alloy bar 

above and below the temperature of transition between the s and n state. In 
particular, we present a description of the model (Section 2), a summary of 
the experimental setup (Section 3), the collected data and analysis (Section 
4) and the comparison between the data and the model (Section 5). 
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2 Discussion of the thermo-acoustic effects 



A pressure pulse is generated in a suspended cylindrical bar in the n state 
following the interactions of an elementary particle with the bulk. This sonic 
pulse, due to the local thermal expansion caused by the warming up, related to 
the energy lost by the particle crossing the material, determines the excitation 
of the vibrational modes of the bar. In the experiment of Ref. [8] an aluminum 
bar was exposed to a proton beam and the theoretical expectations were based 
on a model in which the "amplitude of the fundamental longitudinal mode of 
oscillation" , hereafter referred to as Amplitude, is given by: 

2aLW 



for a beam hitting the center of the cylinder generatrix. In the previous relation 
L, M are respectively length, mass of the cylinder, W is the total energy loss 
of the beam in the bar, a is the linear thermal expansion coefficient and cy is 
the isochoric specific heat. The ratio of the thermophysical quantities a and 
cy is part of the definition of the Griineisen parameter of the material: 

T^^. (2) 

pcv 



where (3 is the volumetric thermal expansion coefficient {(5 — 3q; for alu- 
minum), Kt is the isothermal bulk modulus and p is the mass density. The 
parameter 7 is a very slowly varying function of the temperature when the 
material is in the n state. Solution (1) is a particular case of a more general 
treatment of the problem, which includes the paths of the interacting parti- 
cles in the bulk other than the coordinate of the impact point [9,10,11]. By 
the introduction of a vector field u(x, t) describing the local displacements 
from equilibrium, the amplitude of the mode k of the cylinder oscillation is 
proportional to: 
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where AP*^'^''™ is the pressure pulse due to the sonic source previously de- 
scribed, dW jdx is the specific energy loss of the interacting particle, Jl is 
the cross section of the tubular zone centered on the particle path in which 
the effects are generated and = J dl^V ■ Uk{x)) is a line integral over the 
particle path involving the normal mode of oscillation Uk{x). The Amplitude, 
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as given by (1), can be obtained starting from (3) for a thin bar {R/L <^ 1, 
where R is the bar radius) and for particles hitting the central section. 

In the following, for the material in n state, we will compare the measured 
values of Amplitude to the expected value: 

Xtherm = 5o(l + e) , (4) 



where e is a corrective parameter estimated by a Monte Carlo (MC) simula- 
tion [6], which takes into account the solutions 0[{R/L)'^] for the modes of 
oscillation of a cylinder, the transverse dimension of the beam at the impact 
point and the trajectories of the secondary particles generated in the bar. The 
value of e for the aluminum alloy bar used in the experiment is estimated by 
MC to be -0.04. 

When the material is in the s state, an additional sonic source could be due to 
the local s-n transitions in zones centered around the interacting particle path 
[9,10]. The additional contribution to the amplitude of the cylinder oscillation 
mode k is proportional to: 
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where AV^ and AS are the differences of the volume and entropy in the two 
states of conduction, while A" is the cross section of the tubular zone cen- 
tered on the interacting particle path and switched from s to n state, which 
is given by A" = {dW/dx)/{AT-L/V) [12,13] involving the difference of en- 
thalpy, H., among the two states. The differences can be expressed in terms 
of the thermodynamic critical field and it follows, in first approxima- 
tion, that [14,15]: AV/V={yn-Vs)/V=H^{dHjdP)/{ATi) and AS/V={Sn- 
Ss)/V=-H^{dHjdT)/{A7r). Moreover, by using the difference (A^/\/=(^„- 
Qs) /V=Hl/ {Sti)) of the Gibbs free energy among the two states and by mak- 
ing the hypothesis that He has the parabohc behavior Hc{t) — ifc(0)(l — t^), 
where t = T/T^ and is the transition temperature, it follows that ATi/V — 
H'^{0){1 — t^){l + 3t^)/(87r). In order to compare the observed data with the 
model predictions, we will use the ratio TZ of the contributions to the Am- 
plitude due to local transition effects {Xtrans) and to thermal effects in the n 
state (Xtherm)- Can be expressed as: 
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due to the existent proportionality between the mode amphtude and g. 

In an alternative scenario, which takes into account that local transitions do 
not occur, the Amplitude expected values are given by the relation (4) making 
use of the a and cy values for the s state. 



3 Experimental setup 

The experiment setup has been fully described in Ref . [6] . Here we briefly recall 
that the test mass is a cylindrical bar (R — 0.091 m, L — 0.5 m, M — 34.1 kg) 

made of A15056, the same aluminum alloy (nominal composition 5.2 w% Mg 
and 0.1 w% of both Cr and Mn) used for NAUTILUS. The bar hangs from 
the cryostat top by means of a multi-stage suspension system insuring an 
attenuation on the external mechanical noise of -150 dB in the 1700-6500 
Hz frequency window. The frequency of the fundamental longitudinal mode 
of oscillation of the bar is /o = 5413.6 Hz below T = 4 K. The cryostat 
is equipped with a ^He refrigerator, capable of cooling the bar down to T ~ 
0.5 K. The temperatures are measured inside the cryostat by 10 thermometers 
controlled by a multi-channel resistance bridge. In particular, a calibrated 
Ru02 resistor detects the temperature of one of the bar end faces with an 
accuracy of 0.01 K for T < 4 K. Two piezoelectric ceramics (Pz), electrically 
connected in parallel, are inserted in a slot cut in the position opposite to 
the bar suspension point and arc squeezed when the bar shrinks. In this Pz 
arrangement the strain measured at the bar center is proportional to the 
displacement of the bar end faces. The Pz output is first amplified and then 
sampled at 100 kHz by an ADC embedded in a VME system, hosting the data 
acquisition system. The measurement of the Pz conversion factor A, relating 
voltage to oscillation amplitude, is accomplished according to a procedure 
based on the injection in the Pz of a sinusoidal waveform of known amplitude, 
with frequency /o and time duration less than the decay time of the mechanical 
excitations and on the subsequent measurement of Amplitude. The procedure 
is correct if R/L <^ 1 and a 6% systematic error in the determination of A 
was found. Amplitude is measured according to X = Vq'^^^^ /{GX), where G 
is the amplifier gain and Vf^*^^^ is the maximum of the signal component at 
frequency /o, which is obtained by Fast Fourier Transform algorithms applied 
to the digitized Pz signals. The sign of Amplitude is taken positive or negative 
according to the sign of the first sampling above the noise in the waveform 
generated by the Pz and sampled by the ADC. BTF delivers to the bar single 
pulses of ~ 10 ns duration, containing electrons of 510 ± 2 MeV energy. 
A^e ranges from about 5 x 10^ to 10^ and is measured with an accuracy of 
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~ 3% (for Ne > 5 X 10^) by an integrating current transformer placed close 
to the beam exit point. MC, already introduced in Section 2, estimates an 
average energy lost (AE) ± uae = 195.2 ± 70.6 MeV for a 512 MeV electron 
interacting in the bar and, consequently, the total energy loss per beam pulse 
is given by Vl^ = Ne{AE) , aw = V^Z ctae- 



4 Measurements and data analysis 

Samples of A15056 obtained from the same production batch of the test mass 
have been used to characterize the material at very low temperatures. The 
measurement of the transition temperature to the s state conducted using the 
mutual inductance method gives the value Tc = 0.845 ± 0.002 K and a total 
transition width of about 0.1 K. The smaller value of Tc with respect to pure 
Al (1.18 K) could be ascribed to the presence of Mn impurities in the alloy. In 
fact, experimental studies on AlMn polycrystalline alloys have shown that T^. 
was depressed down to 0.868 K and 0.652 K for Mn concentrations of 440 ppm 
and 900 ppm, respectively [16]. Moreover, A150XX alloys contain inclusions 
of the extremely complex (MgAl) /3 phase [17] and the characterization of su- 
perconducting properties of the alloy /3 — Al3Mg2 shows that — 0.87 K [18]. 
Specific heat data for A15056 are available in literature [19], however, in order 
to completely characterize the production batch, we have performed cy mea- 
surements above and below Tc (Fig. 1) using the calorimetric method of Ref. 
[20]. In the temperature interval 0.9 K < T < 1.5 K the fit of the data points, 
which have an accuracy of 5%, to the function cy/T — T + BT"^ gives the val- 
ues r = 1157 ±31 erg cm^''' K^^ for the electronic specific heat coefficient per 
unit volume in the n state and i? = 0.14 ± 0.01 mJ mol~^ K^** for the lattice 
contribution. If the superconducting properties of A15056 can be described by 
the BCS theory, then Hc{0) ^ 2.42 V^/'^Tc ^ 70 Oe. An independent check of 
the cv{T) behavior is obtained by the measurements at the bar end face of 
the temperature increments due to energy released by each beam pulse. The 
main features of the cy(T) behavior, as obtained by the calorimetry, are well 
reproduced by this method (inset of Fig. (1)). 

The full set of Amplitude measurements (X) normalized to the energy de- 
posited per beam pulse (1^) in the explored temperature interval is shown 
in Fig. 2. For T > 0.9 K, above T^, X has a strict linear dependence on 
W, as expected from the relation (1). The linear fit X = bW (Fig. 3) gives 
b — (2.42 ± 0.17) 10^^° m/J' where the error is determined by the quadrature 
of the beam monitor (3%) and A determination (6%) accuracies. The onset at 
T ~ 0.9 K and the behavior of the superconducting effects are shown in Fig. 2. 
As T decreases, the normalized Amplitude becomes negative, indicating that 
a compression rather than an expansion is generated by the beam interaction 
in the bulk. Its absolute values is greater than b, the normalized Amplitude 
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T[K] 

Fig. 1. A15056 specific heat: calorimetric measurements (5% accuracy). The inset 
shows the calorimetric measurements interpolated by a polynomial (continuous line) 
and the independent cy determination based on the temperature increments at one 
end face of the bar (dots). 



value measured in the n state. The increase in the absolute value of the Am- 
plitude explains the effects seen in cosmic ray observations by NAUTILUS, 
when operated at T = 0.14 K, as due to the conduction state of the material. 

Furthermore, X does not linearly depend on W at fixed T, opposite to what 
has been observed [7] in pure Nb in the s state. The dependence of X/W 
on W in the s state is shown in Fig. (4) representing the data in four non- 
overlapping bands of W. This fact has an impact on the quantification of the 
enhancement of the absolute value of the Amplitude below and far away from 
Tc. Fig. 5 shows the averages of and W in four bins of data collected in 

the temperature interval ranging from 0.55 to 0.60 K, together with the best 
fit given by the exponential {\X/W\) = (8.31 ± 2.88)10"^° g(-26.2±6.3){H/> 
The average energy deposited by the cosmic rays interacting in the NAUTILUS 
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T[K1 

Fig. 2. Measured values of Amplitude [X) normalized to the energy iW) deposited 
in the bar per beam pulse vs. temperature (T) . 

antenna is in the order of 10~^ J [21], much lower than that released in our 
test mass by the beam pulse. We use the ratio = {\X/W\)/b as a factor 
quantifying the Amplitude enhancement in the s state with respect to the n 
one. With reference to Fig. (5), we obtain JF=3.4 ± 1.2 by the extrapolation 
of the fitting exponential to (W) = 10~® J. A value of ~ 3.5 is consistent 
with the cosmic ray observations made by NAUTILUS in the s state [21]. 



5 Comparison with the model 

The Amplitude (X) linearly depends on the deposited energy (W) in the 
model described in Section 2, while a X/W dependence on W is observed in 
the data. Therefore, we try to compare the model predictions to the data in the 
hypothesis that the linear dependence of X on is attained at very low values 
of energy deposition. The application of the model for the expected value 
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0-02 0.04 0,06 0.08 

W[J] 

Fig. 3. T> 0.9 K (n state); Measured values of Amplitude {X) vs. the energy de- 
posited per beam pulse {W). The slope of the fitted line is b = 2.42 10"-'^'^ m/J. 

(Xea;p) computation of the Amplitude in the s state requires the knowledge 
of 1) the thermophysical parameters a„ and cy„ of the material in order to 
evaluate Xtherm for the n state below Tc and 2) the dependence of He on T 
and P for calculating Xtrans via He and its derivatives dHc/ dT and dHd dP. 
The use of relations (1), (2), (4) and (5) allow us to write: 



therm 




i + n) 
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Fig. 4. Data of Fig. (2) ordered in 4 hands of deposited energy. Each hand is identified 
hy (W) ±aw- A dependence of X/W on W can he seen in the s state, e.g. at T=0.6 
K. 



The requirement 1) cannot be fulfilled due to the lack of knowledge of a„ for 
A15056 and we therefore assume that Xtherm/W = b also in the temperature 
interval 0.5 K < T < Tc, due to the fact that 7„, in this interval, is expected to 
have almost the same value as that assumed at slightly higher temperatures. 
In relation to requirement 2), we derive dHc/dT at T < Tc from the He 
parabolic dependence on t, assuming that the unknown dependence of dHc/ dP 
on t at P = for A15056 is equal to that of pure Al. Under this hypothesis, 
dHc/dP can be deduced by interpolating the tabulated values of He as a 
function of T and P contained in Ref. [22]. Inserting numerical values in 
relation (6) gives an average of (Xexp/W) = (— 18±1)10~^° m/J in the interval 
0.55 < T < 0.6 K, where the error does not include systematic contributions 
deriving from the assumptions made. This is to be compared to (X/W) = 
(—8.3 ± 2.8)10~^° m/J, obtained at W = 10~* J from the measurements in 
the same temperature range (see Fig. (5)). This discrepancy can be ascribed 
to the fact that the model, as mentioned by the authors of Ref. [10], considers 
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Fig. 5. 0.55 <T < 0.60 K; averages of the normalized absolute values of Amplitude 
{\X/W\) vs. the the average energy released per beam pulse {W). The line represents 
the fit given by 8.31 10"^° e-26.2(iy> 

superconducting effects in pure materials, wliile tlie intrinsic properties of an 
alloy could determine additional contributions to the expected values of the 
oscillation amplitudes. Indeed, the model has given a satisfactory description 
of the behavior of the experimental data collected with a pure Nb test mass 
in the s state [7]. 

The calculation of X^xp/W for the s state, in the framework of the alternative 
scenario described in Section 2, requires the knowledge of as{= a^ + a'", where 
e and r refer to the electronic and lattice contributions, respectively) and cv,s, 
the former being unknown and the latter measured. The relation (2) gives = 
pcv,s^ll I (^iKt) and = pcy^Y I {'^^t)- Again, for lack of better knowledge, 
we presume the values of the Griineisen parameters and of Kt for AL5056 to 
be similar to those of pure Al. Thus, we use Y = 2.6 [23] in the limit t — 0, 
7| = -11.5 ± 1.0 at t ~ 0.7 [24] and Kt = 79.4 10^ N/m^ near T=0 [25]. The 
insertion of these values in relation (4) gives {Xexp/W) = (—11 ± 1)10^^° m/J 



11 



in the same temperature interval as in the first scenario. As previousy, the error 
is determined only by propagating the errors on the quantities in the right side 
of relation (4). The systematic uncertainties introduced in the calculations of 
the expected values in both scenarios do not allow us to individuate the one 
that better agrees with the experimental data. Moreover, it is interesting to 
derive the predictions at temperature values close to the lower limits of T in the 
H^{P,T) tabulation of Ref. [22]: X^^p/W = (-28±l)10-i° m/J at T = 0.3 K 
in the first scenario and Xg^p/W = (—28 ± 16)10~^° m/J in the alternative 
scenario, where for the latter the large error is due to the uncertainty on 7I 
at this temperature. 

Finally, dissipative efi^ects, which can be inherent in this alloy and due to 
the fiux line motion with consequent entropy transport, could play a role in 
the Amphtude observed values. A clue in this direction hes in the fact that 
Amphtude is not hnearly dependent on W in the s state at fixed T. 



6 Conclusions 

The measurements performed on an A15056 suspended bar, hit by an electron 
beam and operated at temperatures above and below Tc, have shown that in 
the s state, the amplitude of the fundamental mode of the bar is enhanced 
with respect to the n state by a factor ~ 3.5 at T ~ 0.5 K. This factor 
is consistent with the observations made by NAUTILUS on cosmic rays at 
T = 0.14 K. The amplitude change in the s state, following an energetic 
particle interaction, is due to the superconducting properties of the material. 
The absolute value of the normalized Amplitude is enhanced in the AI alloy 
and reduced in Nb. Incomplete knowledge of the involved thermophysical and 
thermodynamic parameters does not allow a full assessment of the model 
describing the underlying physical process in the s state. The effects due to 
cosmic ray interactions could be an important source of noise in future gw 
acoustic and interferometric detectors of improved sensitivities and a complete 
characterization of the thermo-acoustic effects in the test masses operated in 
the s state should be performed by direct measurements of the type shown in 
this letter. 
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